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Abstract 

The action potential signal of nerve and muscle is produced by voltage sensitive channels that include a 

specialized device to sense voltage. Not surprisingly the voltage sensor depends on the movement of 

permanent charges of basic side chains in the changing electric field, as suggested long ago by Hodgkin 

and Huxley and Bezanilla and Armstrong, for example. Gating currents of the voltage sensor are now 

known to depend on the back-and-forth movements of positively charged arginines through the 

hydrophobic plug of a voltage sensor domain. Transient movements of these permanently charged 

arginines, caused by the change of transmembrane potential, further drag the S4 segment and induce 

opening/closing of ion conduction pore by moving the S4-S5 linker. The ion conduction pore is a separate 

device from the voltage sensor, linked (in an unknown way) by the mechanical motion and electric field 

changes of the S4-S5 linker. This moving permanent charge induces capacitive current flow everywhere. 

Everything interacts with everything else in the voltage sensor so everything must interact with everything else 

in its mathematical model, as everything does in the whole protein. A PNP-steric model of arginines and a 

mechanical model for the S4 segment are combined using energy variational methods in which all movements 

of charge and mass satisfy conservation laws of current and mass. The resulting ID continuum model is used 

to compute gating currents under a wide range of conditions, corresponding to experimental situations. 

Conservation laws are partial differential equations in space and time. Chemical-reaction-type models based on 

ordinary differential equations cannot capture such interactions with one set of parameters. Indeed, they may 

inadvertently violate conservation of current. Conservation of current is particularly important since small 

violations (<0.01%) quickly (« 10" 6 seconds) produce forces that destroy molecules. Our model reproduces 

signature properties of gating current: (1) equality of on and off charge in gating current (2) saturating voltage 

dependence in QV curve and (3) many (but not all) details of the shape of gating current as a function of voltage. 

The model computes gating current flowing in the baths produced by arginines moving in the voltage sensor. 

The model also captures the capacitive pile up of ions from the bulk solution in the form of electric double 

layers in vestibules adjacent to both ends of the hydrophobic plug. This pile up is coupled to the movement of 
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arginines, and appears as a spike in the beginning of gating current that have been already observed as the early 
transient gating current in experiments. Our results agree qualitatively well with experiments, and can obviously 
be improved by adding more details of the structure and its correlated movements. The proposed 
continuum model is a promising tool to explore the dynamics and mechanism of the voltage sensor that 
avoids the sampling constraints that so limit the representation of biological concentrations (10 7 - 5xl0 _1 M) 
and currents (time scale > 10" 4 sec) in molecular dynamics simulations that compute atomic motion on the 10" 
15 sec time scale. 
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1. Introduction 

Much of biology depends on the voltage across cell membranes. The voltage across the membrane 
must be sensed before it can be used by proteins. Permanent charges 1 move in the strong electric fields 
within membranes, so carriers of gating and sensing charge were proposed as voltage sensors even 
before membrane proteins were known to span lipid membranes [18-19, 21-23], Movement of 
permanent charges of the voltage sensor is gating current and movement is the voltage sensing 
mechanism. Measurements of gating current were made much easier by biological preparations with 
artificially increased densities of voltage sensors and blocked conduction pores. 

Knowledge of membrane protein structure has allowed us to identify and look at the atoms that 
make up the voltage sensor. Protein structures do not include the membrane potentials and macroscopic 
concentrations that power gating currents, and therefore simulations are needed. Atomic level 
simulations like molecular dynamics do not provide an easy extension from the atomic time scale 2*10" 
16 sec to the biological time scale of gating currents that starts at 50xl0 -6 sec and reaches 50* 10 -3 sec. 
Calculations of gating currents from simulations must average the trajectories (lasting 50* 10" 3 sec 
sampled every 2xl0 16 sec) of ~10 6 atoms all of which interact through the electric field to conserve 
charge and current, while conserving mass. It is difficult to enforce continuity of current flow in 
simulations of atomic dynamics because simulations compute only local behavior while continuity of 
current is global, involving current flow far from the atoms that control the local behavior. It is impossible 
to enforce continuity of current flow in calculations that assume equilibrium (zero flow) under all 
conditions. Periodic boundary conditions are widely used in simulations. Such conditions take a box of 
material and replicate it identically, so the potential at the corresponding edges of the box are identical. 
If the potentials are identical, current will not flow. Periodic boundary conditions of this sort are 
incompatible with current flow from one boundary to the other. Voltage clamp experiments, and natural 
biological function involve current flow from one boundary to another. 

A hybrid approach is needed, starting with the essential knowledge of structure, but computing 


1 Permanent charge is our name for a charge or charge density independent of the local electric field, for example, the charge and charge distribution of Na + but 
not the charge in a highly polarizable anion like Br or the nonuniform charge distribution of H2O in the liquid state with its complex time dependent (and 
perhaps nonlinear) polarization response to the local electric field. 
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only those parts of the structure used by biology to sense voltage. In close packed ('condensed') systems 
like the voltage sensor, or ionic solutions, 'everything interacts with everything else' because electric fields 
are long ranged as well as exceedingly strong [14]. In ionic solutions, ion channels, even in enzyme active 
sites, steric interactions that prevent the overfilling of space in well-defined protein structures are also of 
great importance as they produce long range correlations [17]. Closely packed charged systems are best 
handled mathematically by energy variational methods. Energy variational methods guarantee that all 
variables satisfy all equations (and boundary conditions) at all times and under all conditions and are 
thus always consistent. We have used the energy variational approach developed in [12, 20] to derive a 
consistent model of gating charge movement, based on the basic features of the structure of crystallized 
channels and voltage sensors. The schematic of the model is shown below. The continuum model we use 
simulates the mechanical dynamics in a single voltage sensor, while the experimental data from [7] is 
from many independent voltage sensors. Ensemble averages of independent voltage sensors recording is 
roughly equivalent to macroscopic continuum modeling in a single voltage sensor. It is unlikely one can 
do better until experiments determine the structural source of the variance in the voltage sensor and its 
time, voltage, and chemical dependence on ion type and concentration. 

3. Mathematical model 

The reduced mechanical model for a voltage sensor is shown in Fig. 1 with four arginines Ri, i= 1, 2, 
3, 4, each attached to S4 helix by identical spring with spring constant K. The electric field will drag these 
four arginines since each arginine carries +1 charge. In addition, the charged arginines move the group 
itself. S4 connect to S3 and S5 at its two ends by identical springs with spring constant Kw/2. Once the 
membrane is depolarized from say -90mV inside negative, to +10mV inside positive, arginines together 
with S4 will be driven towards extracellular side. When the membrane is repolarized from say +10mV 
volts to -90mV, the arginines move to the intracellular side. This movement is the basic voltage sensing 
mechanism of ion channel. The movement of S4 triggers the opening or closing of the lower gate mainly 
consisting of S6 of the permeant ion channel by a mechanism widely assumed to be mechanical, although 
we hasten to say that electrical aspects of the linker motion are likely to be involved, as well, perhaps 
acting as a deterministic trigger and control for the otherwise stochastic sudden discontinuous change of 
current seen in single channel measurements of ion flow through the permeation channel. 

When arginines are driven by electric field, they are forced to move through a hydrophobic plug, 
composed by several nonpolar amino acids from SI, S2 and S3 [25], Arginines moving from vestibule, a 
hydrated lumen, on one side, though the hydrophobic plug, and entering into the vestibule on another 
side would involve dehydration when passing through the hydrophobic plug. Therefore, arginines will 
encounter a barrier in the potential of mean force, an energy barrier, mainly dominated by the solvation 
energy [42], Note that Na + and Cl“ (used here as the only solution ions, for simplicity) can only exist in 
vestibules and are not allowed into the hydrophobic plug. The bottoms of the two vestibules on each side 
of the system act as impermeable walls for Na + and Cl“. When the voltage is turned on and off, these two 
walls will form a pair of capacitors storing/releasingnet ion charges in their electric double layers (EDL). 

Molecular dynamics considers every particle individually and must deal with its enormous thermal 
motion (more or less at the speed of sound 1,000 m/s or 1 nm per nanosecond), as it seeks to compute 
gating currents on the biological time scale. We use a continuum treatment that deals with the thermal 
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motion on the average, as temperature or diffusion. Such methods often include the correlations of the 
average electric field, e.g., Poisson-Nernst-Planck (PNP), but not those that differ from atom to atom, as a 
glance at the structure of any protein shows is important. Here we include a steric effect in the PNP model 
[20], and so deal with many of the atomic scale correlations ignored in classical PNP. We hasten to add 
that it is not at all clear how well molecular dynamics deals with correlations produced by the average 
electric field, particularly when they include nonequilibrium components that drive current flow. 
Classical molecular dynamics assumes equilibrium and uses periodic boundary conditions that have 
certain difficulties dealing with the current flow and nonperiodic electric fields found in real devices and 
channels, as we have previously mentioned. Continuum models can adequately describe the time 
behavior of a voltage sensor in action, in which the time span is on the order of several hundreds of 
milliseconds. 

The four arginines Ri, i= 1, 2, 3, 4, are described by their individual density distributions 
(concentrations) a, 7=1, 2, 3, 4, allowing the arginines to interact with Na + and Cl“ in vestibules. The 
density (concentration) distributions represent probability density functions as shown explicitly in the 
theory of stochastic processes, used to derive such equations in [29] using the general methods of [28], 
This practice is used widely in quantum mechanics [3, 34], The important issue here is how well the 
correlations are captured in the continuum model. Some are more likely to be faithfully captured in 
molecular dynamics simulations (e.g., more or less local hard sphere interactions), others in continuum 
models (e.g., correlations induced by far field boundary conditions). 

Here we treat the S4 itself as a rigid body so we can capture the basic mechanism of a voltage sensor 
without considering the full details of structure, which might lead to a three dimensional model hard to 
compute in the time available. We construct an axisymmetric ID model with a three-zone geometric 
configuration illustrated in Fig. 2 following Fig. 1. Zone 1 with z G [0, L R ] is the intracellular vestibule; 
zone 2 with z G [L R , L R + L] is the hydrophobic plug; zone 3 with z G [L R + L, 2L R + L] is the 
extracellular vestibule. Arginines, Na + , and Cl“ can all reside in zone 1 and 3. Zone 2 only allows the 
residence of arginines, albeit with a severe hydrophobic penalty because of their permanent charge, in a 
region of low dielectric coefficient hence called hydrophobic. Spatial variation of dielectric coefficient 
allows construction of many useful devices [15, 41]. Types of boundary conditions and relevant sizes of 
this geometric configuration are also specified in Fig. 2. Note the no-flux boundary conditions specified 
in Fig. 2. One prevents Na + and Cl - from entering the hydrophobic plug (zone 2) with low dielectric 
coefficient. The other boundary condition constrains S4 motion and so prevents the arginines from 
leaving the vestibules into intracellular/extracellular domains. 

We first non-dimensionalize all physical quantities as follows, 


Cl 7 (f> JJ U _ S ~ t =r Dj 9ij 7 Ji 7 7 

Cj =—,&> = —- —, U =-, s = -, t = ——, D; = —, O// =-—, Jj = —-—, / =-, 

1 c 0 ’ ^ k B T/e k B T R R 2 /D x ’ 1 D x ’ ai J k B T/c 0 Jl c 0 D x /r’ ec 0 D x R’ 


where q is concentration of species /, with z'=Na + , Cl , 1, 2, 3, and 4. Each is scaled by c 0 which is the 
bulk concentration of NaCl at the intracellular/extracellular domains. Here c 0 is set to be 184mM, equal 

on both sides, so that the Debye length X D = r r£ ° k ^ T j s lnm when the relative permittivity s r = 80. 

\ c o e 


Note that q, /=!, 2, 3, and 4 needs to satisfy the following additional constraint 


rL + 2LR 

Jo 


A(z)qdz = 1 due 
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to uniqueness of each arginine. <p is the electric potential scaled by k B T/e with k B being the 
Boltzmann constant; T the temperature; e the elementary charge. All relevant external potentials U are 
scaled by k B T. All sizes s are scaled by R, which is the radius of vestibule as shown in Fig. 2. f?=lnm here. 
The time t is scaled by R 2 /D x , with D x being a diffusion coefficient that can be adjusted later to be 
consistent with the time spans of on/off currents measured in experiments (caused by the movement of 
arginines). The diffusion coefficient of species / is scaled by D x . The coupling constant g t j of PNP-steric 
model based on combining rules of Lennard Jones, representing the strength ofsteric interaction between 
species / and j, is scaled by k B T/c 0 [20], Note that here we only consider steric interaction among 
arginines, which is the most significant since arginines are generally crowded in hydrophobic plug and 

vestibules. For simplicity, we assume g t j = j q'^ — j > C7 = 1,2,3,4. The flux density of species i, ] t , is 

scaled by c 0 D x /R, and therefore the electric current / is scaled by ec 0 D x R. For simplicity of notation, we 
will drop ~ for all dimensionless quantities from here on. 

Based on Fig. 2, the governing ID dimensionless PNP-steric equations are expressed below. The first 
one is Poisson equation: 


- L M rA ^) = ^ z ‘ c ‘’ l = Na.au 3,4, 


( 1 ) 


A 2 

with z Na = 1, z ci = —1, Z; = z arg = 1, i = 1,2, 3,4, T = andA(z) beingthe cross-sectional area. For 

zones 1 and 3, T = 1 since here the arginines are fully hydrated with s r = 80. For zone 2, we assume a 
hydrophobic environment with s r = 8, and therefore T = 0.1. 

The second equation is the species transport equation based on conservation law: 

^ + ^(A/i) = 0, i = Na, Cl, 1,2,3,4. (2) 

with the content of/z expressed below based on Nernst-Planck equation for Na + and CP: 

Ji = —Di (j^- + CiZi , i = Na, Cl, z G zones 1 and 3, (3) 

and based on Nernst-Planck equation with steric effect and some imposed potentials for 4 arginines a, 
i= 1, 2, 3 and 4, based on Fig. 2, 

A = _D ‘ (it + to + 'AS + 7?) + >lCl (17 + ^7 + It))- zeallzones, (4) 

^ + f + c 2 (I?+ d ir) + ^ + d ~£ +fr)) ■ z e a “ zones ' & 

h = -°3 (f 1 + z «r e C 3 f + C, + gC 3 ^ + |i)), z € all zones, (6) 

^ = _D -> {it + te + (S + S + gC * + ^ + lr)) ■ z e a11 


zones. 


(7) 


The first and second terms in Eqs. (3-7) describe diffusion and electro-migration term respectively. 
The third term in Eqs. (4-7) are external potential terms with Vi, i= 1, 2, 3 and 4 being the constraint 
potential for the 4 arginines a to S4 represented here by a spring connecting each arginine a to S4 as 
shown in Fig. 1. It is expressed as 
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Vt(z, t) = K(z — (zj + Z S4 (t))) 2 , (8) 

where K is the spring constant, z; is the fixed anchoring position of the spring for each arginine a on S4, 
Z S4 (t ) is the center-of-massz position of S4 by treating S4 as a rigid body. Here we setzi=0.6, Z2=0.2, Z3=- 
0.2, Z4=-0.6 given from the structure that gives the arginine anchoring interval on S4 as 0.4nm. Zs4 (0 
follows the motion of equation based on spring-mass system: 

m s4 dt 2 + b S 4 —^r + K S 4(Z S4 — Z S4 0 ) = Y,t= i K ( z i,cm ~ ( z i + Z s4 )^, (9) 


where m S4> b S4 and K S4 are mass, damping coefficient and restraining spring constant for S4. ■^54,0 is 
the natural position of Zs4 (t) when the net force on the right hand side of Eq. (9) vanishes. Here, z i CM 
is the center of mass for the set of arginines a , which can be calculated by 


_ So + 2 LRA ( z ) zc i dz 
Z i,CM ~ r L + 2LR ., . , > 

/„ a(z)C;dZ 


i=l, 2, 3, 4. 


( 10 ) 


We assume that the spring mass system for S4 is over-damped, which means the inertia term in Eq. (9) 
can be neglected. 

Another external potential in the third term of Eqs. (4-7) is the energy barrier V b , which is nonzero 
only in zone 2 representing mainly the solvation energy barrier from hydrophobicity. Generally, the 
accurate shape of V b in zone 2 requires a calculation of potential mean force (PMF) using molecular 
dynamics simulations by the structure of protein and the possible change in the PMF with membrane 
potential, ionic conditions, consistent with the distribution of charge and boundary conditions, and so on. 
However, here we simply assume a hump shape for the PMF, while we await proper calibrated calculations 
of the PMF and its dependence on membrane potential, etc. 


V b ,max{ tanh(5(z — L R )) — tanh(S(z — L — L R )) — l), z G zone 2, 
0, z G zones 1 and 3. 


( 11 ) 


with V bjnax set to be 5 in the current situation [42], If we set V bmax too large, the gating current would 
be very small since it would be very difficult for arginines to move across this barrier. Note that the tanh 
function is designated to smooth the otherwise top-hat-shape barrier profile, with its awkward infinite 
slopes not good for (numerical) differentiation. It is also generally believed that energy barrier in a 
protein structure does not have a jump. The last term in Eqs. (4-7) is the steric term that accounts for 
steric interaction among arginines [20], Here we set g= 0.5. Larger g implies larger steric effect, but g 
cannot be arbitrarily large due to the limitation of stability. 

Governing equations Eqs. (1-7) were derived by energy variational method based on the following 
energy (in dimensional form): 


E = J ^ cilogci _| V 0| 2 + ^ 


Zj,e Ci<p + ^ 


(V t + V b ) Ci 


arginines 


9ij 

~y~ C i C j 

‘arginines i,j z - 


where the first term is entropy; second and third terms are electrostatic energy; fourth term is constraint 
and barrier potential for arginines; last term is the steric energy term based on Lennard-Jones potential 
[12, 20]. The Poisson equation Eq. (1) is derived by 
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and species flux densities in Eqs. (3-7) are derived by 

SE 

^ ~ 5 ^’ Ji ~ 

where jiq is the chemical potential of species i. 


Pi 

k R T 


qV/q, 


d c ‘ 

Also, here we assume quasi-steady state for Na + and Ch, which means = 0, i = Na, Cl. This is 


justified by the fact that the diffusion coefficients of Na + and Cl - in vestibules are much larger than the 
diffusion coefficient of arginine based on the very narrow time span of the leading spike of gating current 
measured in experiments. The spike comes from the linear capacitive current of vestibule when the 
command potential suddenly rises or drops. This quasi-steady state assumption is essential. Otherwise 
using realistic diffusion coefficients for Na + and Ch would render Eqs. (1-7) too stiff to integrate in time. 
The spike contaminating the gating current can be removed by a simple technique called P/n leak 
subtraction in experiments (see Section 5.4). Howto do this computationally will be discussed in section 
5.3. Boundary and interface conditions for electric potential <p are 

deb deb 

0(0) = v, <P(L~ R ) = 0(z+), T(L- R )A(L- R ) -J- C l- r ) = r(z+)A(z+) -J- (z+), 


0(L r + L-) = ep{h R + L + ), T(L R + L~)A(L r + L-)^(L r + L-) = T(L R + L + )A{L R + L + )^(L R + 

L + ), <p(2L R + L) = 0. (12) 

These are Dirichlet boundary conditions at both ends and continuity of electric potential and 
displacement at the interfaces between zones. Boundary and interface conditions for arginine are 
7i(0, t) = 7i(2L R + L,t) = 0, q(L r , t) = q(L R , t), A(L R )Ji(L R , t) = A(L R )/j(L R , t), 
q(L R +ZT,t) = q(L R +L + ,t), A(L R +L-)J i (L R + L-,t)=A(L R +L + )J i (L R + L + ,t), i = 1,2,3,4. (13) 

No-flux boundary conditions are placed at both ends of the gating pore to prevent arginines and S4 from 
entering intracellular/extracellular domains. The others are continuity of concentration and flux at 
interfaces between zones. Boundary conditions for Na + and Cl - are 
c Wa (0, 0 = c a (0 , t) = c Na (2L R + L,t) = c cl (2L R + L,t) = 1, 

= Jci(k R ,t) = Jnci(L r + L,t) = Jci(L r + L, t) = 0. (14) 

Dirichlet boundary conditions are placed at both ends of the gating pore to describe the concentrations 
for Na + and Cl - are bulk concentration over there and no-flux boundary conditions at both ends of 
hydrophobic plug to describe the impermeability of Na + and Cl“into hydrophobic plug. 

Besides the main input parameter V, which is the voltage bias (command potential) applied, other 
parameters like Di, i= 1, 2, 3,4, K, Ksa, bs4 are also required. Results are especially sensitive to the values of 
K, Ks4, bs4. We have tried and found Di= 50, z=l,2,3,4, K= 3, Ks4=3, bs4=l.S fit best with experiment [7], We 
compare computational outputs with important experiments: (1) gating current versus voltage curve (IV) 
and (2) gating charge versus voltage curve (QV) as well as (3) gating current vs. time curves. 

Usually the electric current in the ion channel is treated simply as the flux of charge and is uniform 
in z when steady. This is not so in the present non-steady dynamic situation, since storing and releasing 
of charge in vestibules and the channel are involved. Here the flux of charge at the middle of hydrophobic 
plug, z= Lr+L/2, was computed to estimate the experimentally observed gating current. However, it is 
actually impossible to experimentally measure the current within the channel since it is impossible to put 
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a probe in the middle of the hydrophobic plug and optical methods to record these variables are so far 
unavailable. In experiments, the voltage clamp technique is used, and on/off gating current through the 
membrane is measured in experiments, which should be equal to the flux of charge at z=0 in the present 
frame work as shown in Fig. 2. The flux of charges at any z position /(z, t) can be related to the flux of 
charges at z=0, 7(0, t), simply by charge conservation: 

y t Qnet(z,t) = 7(0,7) -7(z,t), (15) 

where 

Qnet(z,t ) = (16) 

and flux of charges at any z position 7(z, t) is defined by 

7(z, t) = A (z) Sail iZih{z,t). (17) 

We may as well identify -^Q n et( z <0 as the displacement current, and denote it as 7 disp (z, t), since it will 

be shown later that Eq. (15) is equivalent to Ampere's law in Maxwell's equations, and r^Q n et( z ’ 0 is 

exactly the displacement current in Ampere's law. A general discussion can be found in [10], Hence, Eq. 
(15) can be simply re-written as 

4ot( z < 0 = Kz, t) + Idis V i z ' 0 = 7 (°’ 0, ( 18 ) 

where we define the sum of displacement current and flux of charges as the total current 7 tot (z, t). The 
z-distribution of the total current should be uniform by Kirchhoff s law, which will be computationally 
verified in section 5.3. We are also interested in observing the net charge at vestibules. Taking net charge 
at intracellular vestibule, Q net (.L R , t), as an example, the net charge consists of arginine charges and their 
counter charges formed by the EDL of ionic solution over there. Generally, electro-neutrality is 
approximate but will not be exact there. Note that the existence of EDL's at vestibules, acting as a pair of 
capacitors, is an unavoidable consequence of discontinuities in constitutive properties like permanent 
charge density and dielectric coefficients. Realistic atomic scale simulations will have EDL and capacitive 
currents. More about flux of charge, displacement current and net charge at vestibules will be discussed 
in section 5.3. 

Eq. (15) is consistent with Ampere's law in Maxwell's equations: 

or equivalently, 

V ‘( £ o £ rf+ /) = 0, (20) 

where E is the electric field and J is flux density of charge (current density). Eq. (20) tells us that the 

current is conserved everywhere and it consists of flux of charges J and displacement current s 0 s r —. 

Eq. (20) can be derived by Poisson equation and species transport equation like Eq. (1) and Eq. (2). 
Starting from Poisson equation in dimensional form: 

-V-(£ 0 £ r V$) =p + IjZjeq, (21) 
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or equivalently 


V- (s 0 s r E) = p + 'Z i z i ec i . (22) 

Taking time derivative of Eq. (22), 

v ’( £ “ £ ’-f) = s ' z ' e 7 ? (23) 

and using species transport equation based on mass conservation, 

$ + v-;l = o, (24) 

then 

V ‘ ( £o£r f ) = Zi z i e 5 = “ V ■ Wt = -V ■/. (25) 

which is exactly Eq. (20) by defining 

/ = Ei z i e /i- (26) 

Casting Eq. (20) into the present ID framework and integrate in space from 0 to z, we have 
s 0 £ r A(z) + /(z, t) = £ o £ r A(0) + 1(0, t). (27) 

Comparing with Eq. (18), 

£ 0 £ r A(z) ^ - £ 0 £ r A(0) ^ = / dtop (z, t), (28) 

which justifies the naming of displacement current in Eq. (18). 

To construct the QV curve, we calculate Qi = /„ Lfi A(z) £f =1 qdz, Q 2 = f, Lfi+L A(z) £f =1 qdz, Q 3 = 

f 2L f +L y4(z) Y.i=i c idz, which are the amounts of arginine found in zone 1, 2 and 3, respectively. Usually 

Q 2 ~ 0 due to the energy barrier V b in zone 2. Arginines tend to jump across zone 2 when driven from 
zone 1 to zone 3 when the voltage V is turned on. The number of arginines that move and settles at zone 
3 depends on the magnitude of V. Besides IV and QV curves, the time course of the movement of the 
arginines and of gating charge (and z i CM (t ) and Z S4 (t )), are important to report here. The movement 
of arginines can (almost) be recorded in experiments nowadays by optical methods. Many qualitative 
models accounting for the movement of S4 and conformation change of the voltage sensor have been 
proposed. Readers are referred to the review paper [37] for more details. 

4. Numerical method 

High-order multi-block Chebyshev pseudospectral methods are used here to discretize Eqs. (1-2) 
in space [38], The resultant semi-discrete system is then a set of coupled ordinary differential equations 
in time and algebraic equations (an ODAE system) [2], The ordinary differential equations are chiefly 
from Eq. (2), and algebraic equations are chiefly from Eq. (1) and boundary/interface conditions Eq. (12- 
14). This system is further integrated in time by an ODAE solver (ODE15S in MATLAB [31-32]) together 
with appropriate initial condition. ODE15S is a variable-order-variable-step (VSVO) solver, which is 
highly efficient in time integration because it adjusts the time step and order of integration. High-order 
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pseudospectral methods generally provide excellent spatial accuracy with economically practicable 
resolutions. A combination of these two techniques makes the whole computation very efficient. This is 
particularly important, since numerous computations have to be tried during the tuning of parameters. 


5. Results and discussions 

Here we explored several parameter values to obtain charge movement with kinetics and steady 
state properties similar to the experimentally recorded gating currents. Numerical results based on the 
mathematical model described above were calculated and compared with experimental measurements 
[7], Our one dimensional continuum model has advantages and disadvantages. The lack of three- 
dimensional structural detail means of course that some details of the gating current and charge cannot 
be reproduced. It should be noted however that to reproduce those, one needs more than just static 
structural detail. One must also know how the structures (particularly their permanent and polarization 
charge) change after a command potential is applied, in the ionic conditions of importance. Structural 
detail that does not depend on time can be expected to give an incomplete, probably misleading image of 
mechanism (imagining studying an internal combustion engine without seeing the pistons or valves 
moving). The ID model has the important advantages that it computes the actual experimental results on 
the actual experimental time scale, in realistic ionic solutions and with far field boundary conditions 
actually used in voltage clamp experiments. It also conserves current as we will discuss later 

5.1 QV curve 

When the membrane and voltage sensor is held at a large inside negative potential (e.g., 
hyperpolarized to inside negative say -90mV), S4 is in a down position and all arginines stay in the 
intracellular vestibule. When the potential is made more positive (e.g., depolarized to +10mV, inside 
positive), S4 is in the up position and all arginines are at the extracellular vestibule. 

The voltage dependence of the charge (arginines) transferred from intracellular vestibule to 
extracellular vestibule is characterized as a QV curve in experimental papers and it is sigmoidal in shape 
[7], Fig. 3(a) shows that our computed QV curve— the dependence of Q 3 on V —is in very good 
agreement with the experiment [7], Fig. 3(b) shows the steady-state distributions of Na + , CVand arginines 
in the inside negative, hyperpolarized situation (V=-90mV). As we can see, all arginines stay at 
intracellular vestibule, and none of arginines transferred to extracellular vestibule (Q 3 « 0). 

Fig. 3(c) shows the situation at V=-48mV, which is the midpoint of QV curve. As we can see, each 
vestibule has half of the arginines staying in it (Q 3 = 2): R1 and R2 are expected in extracellular vestibule, 
and R3 and R4 in intracellular vestibule. The concentration distributions of arginines shown in Fig. 3(c), 
interpreted as individual probability density functions, show R1 and R2 residing more in extracellular 
vestibule, and R3 and R4 on the contrary more in intracellular vestibule. Note that there are almost no 
arginines in zone 2 (hydrophobic plug) due to the energy barrier in it. If the command potential at the 
midpoint (Q 3 = 2) is V=0mV, the unforced natural position of S4 is at the middle of gating pore, i.e., 
^S4,0 = Lr + 0.5 L. In experimental measurement [7], Vis actually -48mV instead of OmV for the midpoint. 
This requires Z S40 to be biased from L R + 0.5 L to Z S40 = L R + 0.5L + 1.59lnm. Fig. 3(d) shows the 
situation at full depolarization (V=-8mV). As we can see, all arginines move to extracellular vestibule 
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C<?3 « 4). 


5.2 Gating current 

Fig. 4 shows the time course of gating currents, observed as flux of charge at the middle of gating 
pore /(L r + L/2 , t), due to the movement of arginines when the membrane is largely depolarized, and 
partially depolarized. In the case of large depolarization, V rises from -90mV to -8mV at t=10, holds on 
till t=150, and drops back to -90mV as shown in Fig. 4(a). Time course of gating current and contributions 
of individual arginines are shown in Fig. 4(b). We observe that the rising order of each current component 
follows the moving order of Rl, R2, R3 and R4 when depolarized, and that order is reversed when 
repolarized later. The area under the gating current is the amount of charge moved. Since arginines move 
forward and backward in this depolarization/repolarization scenario, the areas under the ON current 
(arginines moving forward) and the OFF current (arginines moving backward) are same. The areas are 
equal for each component of current as well. The equality of area is an important signature of gating 
current that contrasts markedly with the properties of ionic current. Equality of area has in fact been used 
as a signature that identifies gating current and separates it from ionic current [26, 5]. Since almost all 
arginines move from intracellular to extracellular vestibule when the command potential is suddenly 
made more positive in a large depolarization, the area under each current component is then very close 
to 1. In the case of partial depolarization, V rises from -90mV to -50mV at t=10, holds on till t=150, and 
drops back to -90mV as shown in Fig. 4(c). The time course of gating current and its four components 
contributed by each arginine for this situation is shown in Fig. 4(d). Under this partial polarization, not 
all arginines move past the middle of hydrophobic plug due to weaker driving force in partial 
polarizations compared with large depolarization case. This can be seen in Fig. 4(d), where areas under 
each component current are different, and based on these various areas Rl moves more towards 
extracellular vestibule than R2; R2 more so than R3; and etc. 

The gating currents shown above can be better understood by looking at a sequence of snapshots 
showing the spatial distribution of electric potential, species concentration and electric current. The 
distributions at several times are shown in Fig. 5 for the case of sudden change in command voltage to a 
more positive value, a large depolarization, and Fig. 6 for the case of small positive going change in 
potential, a partial depolarization. Fig. 5 shows that almost all arginines are driven from intracellular to 
extracellular vestibules in case of large depolarization, but only part of the arginines move in case of 
partial depolarization (Fig. 6). The electric potential profiles at t=13 and t=148 show that the profile of 
electric potential as arginines move from left to right even though the voltage is maintained constant 
across the sensor. This is not a constant field system at all! Note that slight bulges in electric potential 
profile exist wherever arginines are dense. This can be easily understood by understanding the effect of 
Eq. (1) on a concave spatial distribution of electric potential. 

In Figs. 5 and 6, total current defined in Eq. (18), though changing with time, is always constant in 
z at all times, which is no surprise because it must satsify Kirchhoffs law. At t=13, a time that gating 
current is not quiescent as shown in Fig. 4(b) and 4(d), we can particularly visualize the z-distributions 
of flux of charges 7(z, t), displacement of current I^spi 2 ’ 0 and total current / tot (z, t) individually in 
Figs. 5 and 6. 
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5.3 Flux of charges at different locations 

Flux of charges 7(z, 7), together with displacement of current I disp (z,t) and total current 
7 tot (z, 7), depicted in Figs. 5 and 6 deserves more discussions here. Flux of charges at the middle of gating 
pore, 7(L R + L/2, 7), and both ends of gating pore, 7(0, 7) and I(2L R + L, 7), should be computed by 

7 (t r + -, 7^ = A (^L r + Yiarginines z iJi (Lr "h ~ (29) 

7(0, 7) = A(0)I 

i=Na,Cl z i Ji (0,7), (30) 

I(2L r + L,t)= A(2L r + L) I i=AffljCl ZjJi(2L R + L, 7). (31) 

Except l(^L R +^,t^j, 7(0,7) and I(2L R + L,t) are trivially zero due to the implement of quasi- 


d c' 

steadiness = 0, i = Na, Cl, in vestibules, which causes J Na and J ci to be uniform in vestibules by Eq. 

(2), and further become zero by the no-flux boundary conditions for Na + and Cl~ at the bottom of 
vestibules as described in Eq. (14). We have to alternatively reconstruct 7(0,7) and 7(2L R + L, t) by 
charge conservation of Na + and Cl~, 

/ (°.0 = A ( z ) Zn 3 ,a z iCidz, (32) 

I{2L r + L, t) = -^j^ R A(z)Y,N a ,a z iCidz. (33) 

After obtaining 7(0, t) and 7(2L R + L, t), we can further reconstruct flux of charges 7(z, t) at zone 1 
(intracellular vestibule) and zone 3 (extracellular vestibule) by charge conservation of Na + and Cl~ 
again, 

7(z,t) = 7(0,7) — ^J 0 Z A(z)X Na ,ciZiQdz, z G [0,L R ], (34) 

7(z,7) = I(2L r +L,t) + yJz Lr+L A( z)Yitia,a z i c idz, z G [L R + L,2L R + L], (35) 

Once flux of charges is analyzed, we can then compute the displacement current based on finding the time 
derivative of Eq. (16). Finally we sum up flux of charges 7(z, 7) and displacement of current 7 d ( Sp (z, 7) 
to obtain total current 7 tot (z, 7) and verify that it satisfies Eq. (18) and is uniform in z. This verification 
is shown in Figs. 5 and 6 at several various times, and is fact true at any time as well. 

In the bottom rows of Figs. 5 and 6 at 7=13, we observe that 7(z, 7) is generally non-uniform in z 
and is accompanied by congestion/decongestion of arginines in between. However 7(z, 7) is almost 
uniform at zone 2 (hydrophobic plug), which means almost no congestion/decongestion of arginines 

occurs there, and therefore no contribution to the displacement current j^Q n eti z > 0 from zone 2. This 

is no surprise since arginines can hardly reside at zone 2 due to energy barrier in it. Because of the energy 
barriers, once arginines leave zone 1 (intracellular vestibule), they immediately jump across zone 2 and 
enter into zone 3 (extracellular vestibule). 

Several things are worth noting in the time courses of 7 (l r + t), and 7(0, 7) illustrated in Fig. 
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7(a) under the case of large depolarization. First, / ( L R + is noticeably larger than 1(0, t ) in ON 
period. This is because their difference, exactly the displacement current 7 d ( Sp , is always negative at zone 
2 when depolarized, since arginines are leaving zone 1 and make ~Qnet < 0 for zone 2. It is expected 

the area under the time course of / ( L R + ^, t) would be very close to 4e, as verified by the time courses 

of Q 3 in Fig. 7(b), while the counterpart area of experimentally measurable 1(0, t ) would be less than 4e 

due to its smaller magnitude compared with / (l R This may account for some experimental 

observations that at most 13e [4, 27, 30], instead of 16e, are moved during full depolarization in 4 voltage 
sensors (for a single ion channel) based on calculating the area under voltage-clamp gating current. 
Therefore, flux of charge at any location of zone 2, though impossible to measure in experiments so far, 
will give us amount of arginines moved during depolarization more reliably than the measurable 1(0, t). 

Second, we see in Fig. 7(a) with magnification in its inset plot that, as in experiments, 1(0, t), but 

not / t has contaminating leading spikes in ON and OFF parts of the current. These spikes are 

linear capacitive currents from solution EDL of vestibules caused by sudden rising and dropping of 
command potential, and are called the early transient gating current in experiments [33, 35-36], In 
voltage-clamp experiments, subtracting this linear capacitive component and removing the spike from 
gating current is done by 'leak subtraction', in various forms, e.g., P/4 (see Section 5.4) In reality, this 
linear capacitive current that is subtracted in this procedure comes from both the lipid bilayer membrane 
in parallel with the channel and vestibule solution in series with the channel. Here, we only considered 
capacitive current from solution EDL of vestibule and ignored the capacitive current of the membrane in 
parallel with the channel (which is actually much larger than vestibule capacitive current) because we 
use Dirichlet boundary conditions for <fi at both ends of gating pore in Eq. (12). Following the idea of 
experiment [7], we calculated 1(0, t) with Arising from -150mV to -140mV at t=10, holding on till t=150, 
and dropping back to -150mV. The reason to choose the rising range of V from -150mV to -140mV is that 
essentially none of the arginines move in this hyperpolarized situation. The voltage step only quickly 
charges and discharges solution EDL in vestibules, and the computational time course of 1(0, t) is just 
two spikes when the command potential rises suddenly and drops suddenly later. Subtracting this 
hyperpolarized 1(0, t), multiplied by a proportion factor, from its original counterpart will then remove 
the spikes, and the unspiked 1(0, t) is shown in Fig. 7(a). 

Third, in Fig. 7(b), the time courses of Qi and Q 3 show that arginines move from the intracellular 
vestibule to the extracellular one when depolarized and move back to intracellular vestibule again when 
repolarized later. As arginines move from one vestibule to another, the concentrations of Na + and Cl" also 
correspondingly change with time at the vestibules, form counter charges through EDL, and balance 
arginine charges at vestibules. However, these movement only maintain approximate, not exact, 
electroneutrality as shown in Fig. 7(b). The violation of electroneutrality is produced by the displacement 
current, that is not negligible we note. 
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As mentioned above, here we used flux of charges at the middle of hydrophobic plug, / ( L R + 
instead of experimentally measurable 1(0, t) to represent‘the gating current' in discussions. This 'gating 
current' just leaves out the displacement current I d i Sp (l R + We use this definition of 'gating 

current' for several reasons: (1) the area under time course of I (^L R +^,t) gives us the amount of 

arginines moved during depolarization more faithfully than 1(0, t) as explained above. The fluxes of 
charge for each arginine shown in Fig. 4 (b) and (d) carry important information about how each arginine 
is moved by the electric field, that will be illustrated in Fig. 8. All these will not be easy to display and 
comprehend if we use 1(0, t) instead. (2) Using 1(0, t) as a definition of gating current would require a 
decontamination by removing the leading spikes in it, which is computationally costly by the procedure 
described above. Especially, it would pose a heavy numerical burden when doing parameter fitting, where 

numerous repeated computations need to be conducted; (3) The shape of / ( L R is actually close 

to the unspiked 1(0, t ). Hence here we used the computationally handy / ( L R + to replace 1(0, t ) 
as the gating current. 

5.4 Comparison with experimental records. 

Our computations have limited fidelity at short times because of time step limitation in integrating 
stiff systems. The spike artifacts are one example, described previously. Experimental measurements [16, 
33] of the fast transient gating current are fascinating but our calculations must be extended to deal with 
them. 

A more general consideration is the subtraction procedure used in experiments to isolate gating 
current from currents arising from other sources. Channels and their voltage sensors are always 
accompanied by large amounts of lipid membranes. Currents through channels and in voltage sensors are 
'in parallel' with large capacitive currents through lipid membranes. The lipid membrane of cells 
introduces a large capacitance ( Cu pid = 6xl0 -7 farads/cm 2 ) that has nothing to do with the capacitive 
(i.e., displacement) currents produced by charge movement in the voltage sensor. Fortunately, this 
capacitance C lipid is a nearly ideal circuit element. Current through the membrane capacitance is entirely 

a displacement current accurately described by i cap = C Upid ^ with a single constant C Upid . V is the 

voltage across the lipid capacitor. Note that i cap does not include any current carried by the translocation 
of ions across the lipid. 

In experimental measurements, this displacement current i cap is always present, in large amounts, 
because voltage sensors (and channels) are always embedded in lipid membranes. Experimental 
measurements mix the displacement current of the lipid membrane and the displacement current of the 
voltage sensor. Indeed, the lipid membrane current dominates the measurement of displacement current 
in native preparations and remains large in systems mutated to have unnaturally large numbers of voltage 
sensors. 
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A procedure to remove the lipid membrane current is needed if the gating current of the voltage 
sensor is to be measured. The procedure introduced by Schneider and Chandler [26] has been used ever 
since in the improved P/4 version developed by Armstrong and Bezanilla [1] reviewed and discussed in 
[8], Also, see another approach in [6, 13]. Schneider and Chandler's procedure [26] estimates the so- 

called linear current i x = C x ^ in conditions in which the voltage sensor and C x behaves as an ideal 

circuit element. An ideal capacitor has a capacitance C x independent of voltage, time, current, or ionic 
composition. The Schneider procedure then subtracts that linear current i x from the total current 
measured in conditions in which the voltage sensor does not behave as an ideal capacitor. The leftover 
estimates the nonlinear properties of the charge movement in the voltage sensor. That is to say it 
estimates the charge movement of the voltage sensor that is not proportional to the size of the voltage 
step used in the measurement. The leftover is called gating current here and in experimental papers. 

The gating current reported in experiments however is missing a component of the displacement 
current of the voltage sensor, because of the Schneider and Chandler's procedure that estimates i x . The 
Schneider and Chandler's procedure does not measure i cap by itself. The linear current i x estimated 
by the Schneider and Chandler's procedure has more in it than the current through the lipid membrane 
capacitor i cap . Rather, it contains the lipid membrane current i cap plus current through any structures 

in the membrane ('in parallel') in which current follows the law i x = C x 

Clearly, some of the current produced by movements of the arginines in the voltage sensor will be a 
linear displacement current, a linear component of gating current. But that component is likely to be 
included in i x and so would not be present in the reported gating current. The linear component of 
gating current is removed by the experimental procedure that is needed to remove the lipid membrane 
displacement current i cap . Other systems may contribute to the linear displacement current as well, for 
example, (1) all sorts of experimental and instrumentation artifacts and (2) displacement current in the 
conduction channel itself. The conduction channel of field effect transistors produces a large capacitance 
described by drift diffusion equations quite similar to the PNP equations of the open conduction channel. 

Our procedure in the computations reported here subtracts a hyperpolarized current with arginines 
not moving in this situation, thereby removing all the currents carried by arginine movement, linear and 
nonlinear in voltage. Most systems have substantial motions that are linear in voltage (even if the system 
is labelled 'nonlinear'). The response of most systems to a voltage step can be written as a sum of terms, 
each with a different dependence on the size of the applied step of voltage V. The first term in that sum 
changes linearly if the voltage step is changed. Higher order terms exhibit a nonlinear dependence on V. 
The linear term is present in most systems, just as it is present in most Taylor expansions of nonlinear 
functions. 

The linear component missed in experiments, and removed in these calculations, may have 
functional and structural significance. The voltage sensor works by sensing voltage, for example, by 
producing a motion of arginines. That motion—the response of the voltage sensor in this model— 
includes a linear component. The signal passed to the conduction channel, to control gating, is likely to 
include the linear component of sensor function, including the linear component of the motion of the 
arginines. Confusion will result if the linear component is ignored when a model is created that links the 
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voltage sensor to the gating process of the conduction channel. 

Direct measurements of the movement of arginines (e.g., with optical methods) are likely to 
include the linear component and so should not agree with experimental measurements of gating current 
or with the currents reported here. 

5.5 Time course of arginine and S4 translocation 

Time course of Q (amount of arginines moved to extracellular vestibule, equal to Q 3 here) and 
center-of-mass trajectories of individual arginines [z iCM , i= 1, 2, 3, 4) and S4 segment (Z S4 ) are shown in 
Fig. 8 with (a) and (b) for the case of large depolarization and (c) and (d) for the case of partial 
depolarization. In the case of large depolarization, from Fig. 8(b), we can see arginines and S4 are driven 
towards extracellular vestibule once the electric field is turned on. Their z-positions quickly reach 
individual steady states with almost all arginines transferred to extracellular vestibule as previously 
shown in Fig. 5 and therefore Q is close to 4 as shown in Fig. 8(a). Arginines and S4 move back to 
intracellular vestibule once the voltage drops back to -90mV. From Fig. 8(b), the forward moving order of 
arginines is Rl, R2, R3 and R4, and the backward moving order is on the opposite R4, R3, R2 and Rl. This 
obviously coincides with the structure. Note that S4 is initially farthest to the right but lags behind Rl and 
R2 during movement in depolarization as shown in Fig. 8(b). This is certainly because S4 is relaxed to an 
almost unforced situation close to its natural position ■^54,0 ■ We can further calculate the displacements 
of each arginine and S4 during this full depolarization, and find A z 1CM ~ Az 2i cm ~ A z 3 CM ~ 1.93nm, 
Az 4 CM =1.76nm, AZ 54 =l,51nm. Besides almost the same displacements for Rl, R2 and R3, their average 
moving velocities are also very close to each other. This seems to suggest a synchronized movement 
among Rl, R2 and R3. Also, we can see the movements of arginines contribute significantly to the 
movement of S4 segment. This can be seen by steady state z-position of S4 derived from Eq. (9), 
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Experimental estimates of S4 displacement duringfull depolarization ranged from 2-20 A [24, 37], 
This large deviation accommodates several models for voltage sensing: transporter model, helical screw 
model and paddle model [37], Our AZ 54 =1.51nm here is large enough that seems to agree better with 
experimental estimates requiring large displacements, such as the paddle model. However, our model is 
more based on helical screw model, which is known to have shorter displacements. A plausible 
explanation for our over-estimated A Z S4 is that our ID model has made the known helical path into a 
linear path. It would be shorter for displacement perpendicular to the plane of the membrane. The 
diameter of an alpha helix between alpha carbons is approximately 1 nm, therefore if the rotation were 
90°, the displacement in ID would be extended by 0.5x0.5*71=0.78 nm, making the linear translation 
perpendicular to the membrane only 0.73 nm. This number is very close to the displacement 
perpendicular to the membrane that is estimated when comparing the open-relaxed state crystal 
structure of Kvl.2 [9] and the consensus closed structure that has been derived from experimental 
measurements [40], 

For partial depolarization, arginines and S4 are the same driven towards extracellular vestibule 
once the electric field is turned on and return when voltage drops. However, since the driving force is 
weaker than in a large depolarization, their z-positions do not reach steady states as they do during a full 
depolarization before they return to intracellular vestibule once the voltage drops. This behavior can be 
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seen from Fig. 8(c), that Q reaches 1.57 at most which should be 2 instead if equilibrium, which can be 
reached if time is long enough, as shown in QV curve of Fig. 3(a). Fig. 8(d) shows that S4 initially farthest 
to the right lags behind R1 during movement and is almost caught up by R2. The maximum displacements 
of arginines and S4 calculated from Fig. 8(d) are A z 1CM = 1.36nm, A z 2CM = 0.966nm, 
Az 3 cm =0.459nm, Az 4CM =0.316nm, and A Z 4 CM =0.616nm. 

5.6 Family of gating currents for a range of voltages 

Fig. 9(a) shows the time courses of gating current for a range of voltages V, ranging from -62mV to 
-8mV. We can see the area under gating current, for both on and off parts, increases with V since more 
arginines are transferred to extracellular vestibule as V increases. However, this area will saturate with 
further increasing V since all arginines would then be transferred to the extracellular vestibule. The 
shapes of this family of gating currents agree well with experiment [7] in both magnitude and time span. 
We can characterize the time span by fitting the time course (only the decay part) of a gating current by 
ae~ t/fTl + be -t / T2 , x 1 < t 2 , as generally conducted in experiment [7], where t 1 is the fast time constant 
and t 2 is the slow time constant. Usually the movement of arginines is dominated by t 2 . Here t 2 was 
calculated from simulation and compared with experiment [7] as shown in Fig. 9(b). Since in our 
computation the time is in arbitrary unit, we have scaled the time to have the maximum t 2 to fit with its 
counterpart in experiment [7], The trend of t 2 versus V in our result agrees well with experiment [7], 
For the left branches to the maximum point in Fig. 9(b), simulation result fits very well with experiment. 
For the right branches to the maximum point, simulation result overestimates t 2 compared with 
experiment. This is consistent with the observation that the amount of transferred charges Q saturates 
faster in experiment data than in present simulation as ^increases as shown in the QV curve of Fig. 3(a). 
This phenomenon is related to cooperativity of movement among arginines, that will be further discussed 
below. 

5.7 Effect of voltage pulse duration 

Fig. 10 shows the effect of voltage pulse duration with Fig. 10(a) for the case of partial 
depolarization and Fig. 10(b) for the case of full depolarization. Magnitude and time span of gating 
current are influenced by pulse duration in both cases, but the shape will asymptotically approach the 
same as pulse duration increases. This behavior occurs because driving arginines towards extracellular 
vestibule by the command potential takes time. If the pulse duration is long enough, the time course of Q 
will approach its steady state like in Fig. 8(a) for large depolarization. Partial depolarization takes longer 
time to reach its steady state as demonstrated in Fig. 8(c). Again these shapes of gating currents in Fig. 10 
compare favorably with experiment [7], 

6. Conclusions 

The present ID mechanical model of the voltage sensor tries to capture the essential structural 
details that are necessary to reproduce the basic features of experimentally recorded gating currents. 
After finding appropriate parameters, we find that the general kinetic and steady-state properties are 
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well represented by the simulations. The continuum approach seems to be a good model of voltage 
sensors, provided that it takes into account all interactions, and satisfies conservation of current. The 
continuum approach used here describes the mechanical behavior of a single voltage sensor, but 
experimental measurements are ensemble averages of multiple ion channels (and hence multiple voltage 
sensors). There is no conflict, however. Experiments average of measurements of large numbers of 
voltage sensors as does our macroscopic model. 

Here we simplify the profile of energy barrier in hydrophobic plug since the potential of mean 
force (PMF) in that region, and its variation with potential and conditions, is unknown. Therefore the next 
step is to model the details of interactions or the moving arginines with the wall of the hydrophobic plug. 
There is plenty of detailed information on the amino acid side chains in the plug and how each one of 
them has important effects in the kinetics and steady-state properties of gating charge movement [25], 
The studied side chains reveal steric as well as dielectric interactions with the arginines that the present 
model does not have. On the other hand, the power of the present mathematical modeling is precisely the 
implementation of interactions, therefore we believe that when we add the dielectric details of the 
channel a better prediction of the currents should be attained. Also, here the plug energy barrier 14 is 
assumed to be independent of time. However, once the first arginine enters the hydrophobic plug by 
carrying some water with it, this partial wetting of pore will lower 14, chiefly consisting of solvation 
energy, and enable the next arginine to enter the plug with less difficulty. This might explain the 
cooperativity of movement among arginines when they jump through the plug. This remains to be 
included in the future modeling. Besides, it has been reported that a very strong electric field might affect 
the hydration equilibrium of the hydrophobic plug and would lower its hydration energy barrier as well 
[39], This cooperativity of movement may help explain the quick saturation in the upper right branch of 
QV curve (and smaller T 2 ). This also remains to be considered in the future modeling. 

Also, further work must address the mechanism of coupling between the voltage sensor 
movements and the conduction pore. It seems likely that the classical mechanical models will need to be 
extended to include coupling through the electrical field. It is possible that the voltage sensor modifies 
the stability of conduction current by triggering sudden transitions from closed to open state, in a 
controlled process reminiscent of Coulomb blockade. 
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Figure 1. Geometric configuration of the model including the attachments of arginines to the S4 segment. 



Figure 2. Following Figure 1, an axisymmetric 3-zone domain shape is designated in r-z coordinate for 
the current ID model. Here the diameter of hydrophobic plug is 0.3nm (arginine's diameter); L=0.7nm; 
Lff=1.5nm; radius of vestibule isff=lnm. BC means Boundary Condition 
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Figure 3. (a) QV curve and comparison with [7], Steady-state distributions for Na, Cl and arginines at (b) 
k=-90mV, (c) k=-48mV, (d) k=-8mV. 
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Figure 4. (a) Time course of V rising from -90mV to -8mV at t=10, holds on till t=150, and drops back to - 
90mV, (b) time course of gating current and its components corresponding to (a), (c) time course of V 
rising from -90mV to -50mV at t=10, holds on till t=150, and drops back to -90mV, (d) time course of 
gating current and its components corresponding to (c). 
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Figure 5. The four panels on the top row are species concentration distributions at t=0, 13, 148, for the 
case of large depolarization with Arising from -90mV to -8mV at t=10, holding on till t=150, and dropping 
back to -90mV. The four panels on the middle row are concurrent electric potential profiles. The four 
panels on the bottom row are concurrent electric current profiles with components of flux of charge, 
displacement current and total current. 
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Figure 6. The four panels on the top row are species concentration distributions at t= 0, 13, 148, for the 
case of large depolarization with V rising from -90mV to -50mV at t=10, holding on till t=150, and 
dropping back to -90mV. The four panels on the middle row are concurrent electric potential profiles. The 
four panels on the bottom row are concurrent electric current profiles with components of flux of charge, 
displacement current and total current. 
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Figure 7. (a) Time courses 
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and unspiked 1 (0, t) for the case of full 
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depolarization with V rising from -90mV to -8mV at t= 10, holding on till t=150, and dropping back to - 
90mV. The inset plot is a blow-up of ON-current to visualize the difference of 1(0,t) andunspiked 1(0,t) 

more clearly, (b) Time courses of Q t , Q 3 , Jq R (c Ni1 — c cl ) dz, and J^ LR ^ L (c Na — c cl ) dz under the same 


depolarization scenario as (a). 



Figure 8. (a) and (c) are time courses of amount of arginines moved to extracellular vestibule, (b) and 
(d) are center-of-mass trajectories of individual arginines and S4. (a) and (b) are the case of large 
depolarization with V rising from -90mV to -8mV at t=10, holding on till t=150, and dropping back to - 
90mV. (c) and (d) are the case of partial depolarization with V rising from -90mV to -50mV at t=10, 
holding on till t=150, and dropping back to -90mV. 



Figure 9. (a) Time courses of gating current with voltage rising from -90mV to FmV at t= 10, holds on till 
t=150, and drops back to -90mV, where V=-62, -50, ... -8 mV. (b) t 2 versus V compared with experiment 
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Figure 10. Effect of voltage pulse duration: (a) V increases from -90mV to -35mV at t=10 and drops back 
to -90mV at various times, (b) V increases from -90mV to OmV at t=10 and drops back to -90mV at 
various times. 


26 



















